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CHAPTER  1 


INTRODUCTION 


1 . 1  Objective  and  Scope 

The  availability,  diversity,  and  utilization  of  anisotropic  materials 
have  increased  at  a  remarkable  rate  during  the  past  few  decades.  When 
attempts  were  made  to  predict  the  behavior  of  these  materials,  It  was 
evident  that  the  available  theoretical  methods  were  far  from  satisfactory. 
As  a  result,  much  renewed  Interest  has  been  directed  to  the  theory  of  anis~ 
troplc  elasticity.  In  particular,  the  problem  of  stress  concentration 
around  cavities  In  a  medium  occupying  the  entire  space  has  long  been  the 
subject  of  numerous  investigations,  not  only  because  of  its  important  role 
in  fundamental  elasticity  problems  but  also  because  of  modern  engineering 
design  concerns. 

The  objective  of  this  study  is  to  develop  an  explicit  analytical  solu¬ 
tion  for  the  non-axlsymraetric  problem  of  three-dimensional  stress 
concentration  in  a  transversely  isotropic  medium  containing  a  spheroidal  or 
spherical  cavity. 

1 . 2  History  and  Literature  Survey 

Stress  concentration  can  be  described  as  the  local  intensification  of 
stress  as  the  result  of  nearby  changes  in  geometry  or  discontinuities  in 
material  or  load.  When  a  cavity  exists  inside  the  material,  it  gives  rise 
to  stress  concentration  and  often  leads  to  structural  failure.  Since  the 


first  investigations  of  the  stress  concentration  problem  during  the  second 
half  of  the  19th  century,  considerable  research  has  been  done  on  the 
problems  of  determining  the  stress  concentration  produced  in  an  elastic 
body  by  a  spherical  or  spheroidal  cavity.  A  general  review  of  stress 
concentration  in  linear  elasticity  has  been  given  in  a  1958  survey  paper  by 
Sternberg  [40].  Eight  years  later,  Neuber  and  Hahn  [31]  again  reviewed 
stress  concentration  in  scientific  research  and  engineering.  Because  of 
the  enormous  number  of  articles  concerning  this  subject,  we  will  only 
briefly  present  the  work  relevant  to  the  three-dimensional  stress  concen¬ 
tration  problems  in  an  infinite  medium  containing  a  spheroidal  or  spherical 
cavity  or  inclusion. 

As  early  as  1866,  Lame  [22]  investigated  the  problem  of  the  spherical 
shell  under  uniform  external  pressure.  A  few  decades  later,  Larmor  [24] 
studied  the  effect  of  a  spherical  cavity  on  a  field  of  pure  shear.  A 
closed  form  solution  for  a  spherical  cavity  in  an  infinite  isotropic  medium 
under  uniaxial  tension  was  presented,  without  derivation,  by  Southwell  and 
Gough  [39].  Goodier  [14]  investigated  the  concentration  of  stress  around 
spherical  and  cylindrical  inclusions  and  obtained  in  a  general  form  the 
Southwell  and  Gough  solution.  The  spheroidal  cavity  under  uniform  axial 
tension,  pure  shear,  and  torsion  was  first  considered  by  Neuber  [30].  An 
investigation  of  the  stress  concentration  around  an  ellipsoidal  cavity  in 
an  isotropic  medium  under  arbitrary  plane  stress  perpendicular  to  the  axis 
of  revolution  of  the  cavity  was  undertaken  in  1947  by  Sadowsky  and 
Sternberg  [35].  They  further  presented  an  analysis  of  the  triaxial  ellip¬ 
soidal  cavity  in  an  infinite  medium  under  a  uniform  stress  field  at 
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Infinity  [36],  In  1951,  Edwards  [9]  obtained  a  solution  for  the  stress 
concentration  around  a  spheroidal  inclusion  in  a  medium  subjected  to  a  uni¬ 
form  stress  field  at  infinity.  A  general  representation  of  the  stress 
concentration  problems  in  isotropic  materials  appeared  in  1965  when 
Podil'chuk  [32]  studied  the  deformation  of  an  axi symmetrically  loaded  elas¬ 
tic  spheroidal  cavity.  Later,  he  extended  his  study  to  include  the  non- 
axisymmetric  deformation  of  spheroidal  cavities  [331. 

Although  the  stress  concentration  problem  has  been  extended  to  aniso¬ 
tropic  materials,  very  few  three-dimensional  anisotropic  elasticity 
problems  have  been  solved.  Chen  [51  solved  the  problem  of  uniaxial  axisym- 
metric  tension  applied  at  large  distances  from  a  cavity  in  a  transversely 
isotropic  medium.  In  1971,  he  also  presented  the  general  solution  for  an 
infinite  elastic  transversely  isotropic  medium  containing  a  spheroidal 
inclusion  [7]  with  the  restriction  that  the  prescribed  stress  field  is 
axisymmetric  and  torsionless.  In  addition,  he  solved  the  spheroidal  inclu¬ 
sion  problem  under  pure  shear  in  and  out  of  the  plane  of  isotropy  and 
modified  Bose's  solution  [3]  for  the  torsion  of  a  transversely  isotropic 
medium  containing  an  isotropic  inclusion,  to  the  case  of  a  spheroidal 
inclusion  in  which  the  material  may  be  transversely  isotropic,  provided 
that  the  axes  of  anisotropy  of  both  the  medium  and  the  inclusion  be  the 
same  (4).  To  the  the  best  of  our  knowledge,  Chen's  solution  for  the  stress 
concentration  around  spheroidal  inclusions  and  cavities  in  a  transversely 
isotropic  material  under  pure  shear  are  the  only  non-axy isymmetr ic  problems 
whose  solutions  are  available.  The  present  work  extends  the  work  by  Chen 
[4-7 |  whose  solutions  can  be  deduced  from  this  approach. 


-v  .v  a‘ ■  *.  v.'  o  .  **  i  . 
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1.3  Organization  of  the  Study 

After  the  aforementioned  historical  review  of  the  problem  of  stress 
concentration,  the  basic  formulae  and  the  potential  functions  approach  in 
the  three-dimensional  elasticity  theory  of  transversely  isotropic  materials 
are  presented  in  Chapters  2  and  3,  respectively. 

In  Chapter  4,  we  introduce  the  special  coordinate  systems  used  to 
solve  the  cavity  problem  and  the  potential  functions  in  terms  of  the  asso¬ 
ciated  Legendre  functions  of  the  first  and  second  kind  as  well  as  some 
important  identities  that  the  potential  functions  satisfy.  The  formulation 
and  the  analytical  approach  of  solving  the  first  and  second  boundary  value 
problems  of  elasticity  theory  for  a  spheroidal  cavity  embedded  in  a  trans¬ 
versely  isotropic  medium  are  presented  in  Chapters  5  and  6,  Problems  for 
cavities  and  inclusions  subjected  to  a  general  constant  stress  field  at 
infinity  are  solved  at  the  end  of  these  two  chapters. 

In  Chapter  7,  numerical  investigation  of  the  stress  concentration 
factors  associated  with  axisymmetric  and  non-axisymmetric  problems  for  a 
variety  of  materials  is  carried  out.  The  effect  of  anisotropy  on  the 
stress  concentration  factor  is  discussed. 

Chapter  8  summarizes  the  developments  of  this  study  and  makes  recom¬ 


mendations  for  further  study. 


CHAPTER  2 


TRANSVERSELY  ISOTROPIC  MATERIALS 


2.1  Definition 


A  transversely  isotropic  material  is  an  anisotropic  one  for  which  the 
Hookean  matrix  at  a  point  remains  invariant  under  an  arbitrary  rotation 
about  an  axis  (the  "axis  of  elastic  symmetry"  of  the  material). 

There  are  some  crystalline  materials  recognized  as  transversely 
isotropic.  These  Include  ice,  cadmium,  cobalt,  magnesium,  and  zinc.  Other 
examples  of  transverse  isotropy  are  provided  by  materials  having  hexagonal 
structures  such  as  fiber-reinforced  composites  [2]  as  shown  in  Fig.  1  and 
the  honeycomb  structures  of  Fig.  2,  as  well  as  laminated  media  (Fig.  3). 
Stratified  rocks  and  soils  can  also  be  modeled  as  a  homogeneous  transverse¬ 
ly  Isotropic  medium  [23,37,38,42]. 

2.2  Fundamental  Formulae 

Consider  a  homogeneous  transversely  isotropic  elastic  medium  occupying 
a  region  of  the  three-dimensional  Euclidian  space  referred  to  a  fixed 
cylindrical  coordinates  system  (r,  9  ,  z)  in  which  the  z-axis  coincides  with 
material  axis  of  symmetry.  Let  (ur,  uQ ,  uz )  denote  the  cylindrical  scalar 
components  of  the  displacement  vector  and  (err,  ee0  ,  yzz,  yqz  ,  Yrz,  Yrg) 
denote  engineering  strains  by 
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In  the  absence  of  electrical  and  thermal  effects,  the  generalized  Hooke’s 
law  [16,26)  can  be  expressed  In  the  form: 


r 

u 

u 

D 

1 

a00 

°ZZ 

°0Z 

°rz 

Q 

CD 

1 _ 

11 

C12 

C13 

0 

0 

0 

e 

!  " 

12 

C11 

C13 

0 

0 

0 

! «-  ! 

1  -  -  | 

i  i 

13 

C13 

C33 

0 

0 

0 

i 

ezz  i 

1 

0 

0 

0 

c44 

0 

0 

i 

V  1 

-z  1 

1 

0 

0 

0 

0 

C44  ° 

i  i 

!  ’  rz  ! 
i  i 

0 

0 

0 

0 

0 

2  (C11  "  C12) 

Lv ; 

c13  ’ 

c33 

and  c44 

are 

the 

five  Independent  elastic 

1  (2.2.2) 


characterizing  the  medium  and  (orr,  °ee»  °zz»  °nz »  °rz»  °re)  are  the 
components  of  the  stress  tensor, 
c, 
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We  may  also  write  for  convenience 


In  the  literature,  tba  coefficients  a^,  a^»  *jj  »  *33*  and  *44  are  called 

"compliances"  vharaaa  the  elaatlc  conatanta  are  called  "self fnesses".  We 

nay  alao  write,  for  convanlanca,  a.,  ■  2(a, ,  -  a,,). 

00  XI  12 

It  la  evident  from  Eqa.  (2.2.2)  and  (2.2.3)  that  the  elastic  constants 
(atlffnaases)  can  be  converted  to  the  compliances  and  vice  versa  by  the 
standard  determinant  procedure  for  solving  simultaneous  equations.  By 


doing  so,  the  following  relationships  are  found: 
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2 . 3  The  Engineering  Constants 


Soae  authors  replace  Che  elastic  constants  and  the  compliances  by  the 
so-called  engineering  constants  which  can  be  Interpreted  to  be  Young's 
moduli,  the  shear  moduli,  and  Poisson's  ratios  associated  with  various 
directions  (see  [25)).  If  we  Introduce  these  constants  which  are  related 
to  the  compliances  by: 
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Equations  (2.2.3)  cake  the  form: 


www>viwvrt,»’. 
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where: 

E  and  E  are  Young's  moduli  In  Che  plane  of  isocropy  and  perpendicular 
co  lc  respectively. 

'•>  is  Poisson's  racio  characterizing  transverse  contraction  in  the 

plane  of  isocropy  when  tension  is  applied  in  chis  plane. 

•j  is  Poisson's  ratio  characterizing  transverse  contraction  in  the 

plane  of  Isocropy  when  tension  is  applied  in  a  direction  normal  co  the 
plane  of  Isocropy. 

G  and  G  are  shear  moduli  for  Che  plane  of  isocropy  and  any  plane  per¬ 
pendicular  co  lc  respeccively. 


the  engineering  conscancs  and  Che  elasdc  constants  can  be  related  co 
each  ocher  by  means  of  Eqs.  (2.2.5)  and  (2.3.1).  These  relations  are  given 
by: 
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For  lsocroplc  aaterlals,  the  engineering  constants  and  the  elastic  con¬ 
stants  are  reduced  to 
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Inequalities 

ratio: 


(2.4.3)  become  for  materials  having  posiclve  Poisson's 


E>0  ,  E  >  0  ,  C  >  0  ,  0  <  v  <  1  -  2"2  - 

E 


the  last  inequality  of  (2.4.4)  is  plotted  in  Fig.  4.  The  interior  area 

E 

bounded  by  a  parabola  corresponding  to  a  specific  value  of  -  ,  represents 

h 

the  valid  domain  of  variation  of  Poisson's  ratio.  Line  01  represents  the 

E 

isotropic  materials.  It  is  noteworthy  that  for  values  of  ^  >  2,  Poisson's 
ratio  v  might  exceed  unity. 

2.5  Values  of  the  Elastic  Constants 

The  numerical  values  of  the  elastic  constants  for  a  variety  of  differ¬ 
ent  hexagonal  crystalline  and  non-crystalline  materials  are  reported  by 
Hearmon  in  the  revised  edition  of  Landolt-Bornstein  [19].  Some  of  these 
materials  are  listed  in  Tables  1  and  2.  In  both  tables,  the  unit  used  for 
the  elastic  constants  c^j  is  the  gigapascal  (GPa  "  10^  dyne/cm^)  and  the 
unit  used  for  the  compliances  a^j  is  the  reciprocal  of  terapascal 
(TPa  -  103  GPa). 

Moreover,  the  elastic  constants  of  fiber  reinforced  materials  can  be 
calculated  by  utilizing  the  technique  developed  by  Hlavacek  [17,18]  or  by 
Hashin  and  Rosen  [15].  Achenbach  [1]  has  presented  a  procedure  for  the 


computation  of  the  effective  moduli  of  laminated  media  while  Salaraon  [37] 
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derived  expressions  for  the  elastic  moduli  of  a  stratified  rock,  mass  and 
Wardle  and  Gerrard  [42]  discussed  the  restrictions  on  the  ranges  of  some 
of  the  five  elastic  constants. 


V 


Zirconium,  a-Zr  144  74  67  166  33.4  10.1  -4.0  — 2.4  8.0 


Table  2.  Elastic  Constants  for  Hexagonal  Systems 
Non-Crystalline  Materials 
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CHAPTER  3 


POTENTIAL  SOLUTIONS  FOR  TRANSVERSELY  ISOTROPIC  MATERIALS 

3.1  Introduct ion 

Ac  Che  beginning  of  Che  twentieth  century,  Michell  [28)  and  Fredholm 
[13)  generaCed  important  three-dimensional  solutions  in  transversely 
isotropic  materials.  Lekhnitskii  (25)  presented  a  solution  procedure  for 
ax i symmetric  torsionless  problems  in  terms  of  a  single  stress  function 
satisfying  a  fourth-order  partial  differential  equation.  Elliott  [10-11) 
presented  a  comparable  solution  in  terms  of  two  stress  functions  each 
satisfying  a  second-order  partial  differential  equation.  Shortly  after 
Elliot's  presentation,  Eubanks  and  Sternberg  [12)  proved  the  completeness 
of  the  Lekhnitskii  and  Elliot  representations  for  the  case  of  axisymmetry. 
Hu  [21]  introduced  a  third  potential  function  suitable  for  non-axisymmetric 
stress  field.  This  solution  was  implied  in  Fredholm's  work  [13).  The 
completeness  of  the  three-function  approach  for  general  non-axisymmetric 
problems  was  proven  in  Reference  [41], 

3.2  Fundamental  Formulae 

Adopting  the  notations  in  Section  2.2  and  combining  Eqs.  (2.2.1)  and 
(2.2.2),  we  obtain  the  stress-displacement  relations  that  take  the  form 
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It  has  been  shown  [41 j  that  In  the  absence  of  body  forces,  the  following 
representation  is  complete 
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where  <t>^(r,9,z),  and  i>(r,9,z)  satisfy  the  equations 
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In  ch«s«  equations,  k^  sad  k^  art  the  roots  of  che  equation 

c44(c13+C44)k2+[(c13+C44)2+c44‘CllC33lk+C44(c44+cl3)  *°  (3-2‘5) 
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and  and  are  che  roots  of  Che  equation 
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wheraas  is  defined  by 
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(3.2.7) 
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The  conscancs  and  are  either  real  or  complex  conjugate  (with  a  real 

pare  different  from  zero)  depending  upon  the  elastic  constants,  but  che 
2 

constant  1«  always  real  and  positive.  We  also  specify  chat  v^,  v^,  and 
always  have  positive  real  parts. 

2  2 

The  constants  k^  and  kj  are  related  to  and  respectively  by 
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or,  equivalently 
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By  using  Eqs.  (3.2.5),  (3.2.6)  and  (3.2.8),  1c  is  noc  difficult  to  obtain 
the  following  identities 
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Sometimes,  it  is  convenient  to  employ  three  new  variables  defined  by 
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These  new  three  variables  generates  three  distinct  spaces  (r,  ■:,*  ), 

(r,9,z  )  and  (r,3,z^)  different  from  the  physical  space.  The  differential 
2 

operations  7^  defined  by  Eq.  (3.2.4)  become,  in  the  spaces  (r, •),*.),  the 
Laplaclan  defined  by 


(3.2.12) 


the  displacement  components  in  Eq.  (3.2.2)  can  be,  alternately,  expressed 
in  the  following  form 
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and,  consequently,  the  stress  components  can  be  written  In  the  form 
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Table  4.  Values  of  \> .  and  k.  for  Selected  Materials* 

C>  1 _ ] _ 

■v  2  2 

Material  v,  v,  k 
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Berylium 

1.216 

1.050 
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0.3058 

1 

0.8491 

+ 

0.5283 

1 

Cadmium 

0.5297 

0.575 

+ 

0.3283 

i 

0. 7’79 

+ 

0.6284 
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Zinc 

0.5915 

0.286 

+ 

0.5411 

1 

0.0847 

0.9964 

i 

CHAPTER  4 


THE  SPHEROIDAL  CAVITY  PROBLEM 


4.1  Statement  of  Problem 


Consider  an  infinite  elastic  transversely  isotropic  medium,  as  shown 
in  Fig.  5,  containing  a  spheroidal  cavity  whose  surface  is  defined  by 


where  a  and  b  are  the  two  semi-axes  of  the  spheroid. 

In  the  absence  of  body  forces,  w#  wish  to  develop  displacement  fields 
whose  scresses  vanish  aC  infinity  and  which  satisfy  given  boundary  condi¬ 
tions  on  the  cavity  surface.  In  this  chapter,  we  introduce  applicable 
coordinate  systems  and  potential  functions.  General  solutions  to  problems 
in  which  either  the  displacement  vector  or  the  surface  traction  vector  is 
prescribed  on  the  surface  of  the  cavity  will  be  considered  in  the  following 
chapters. 

^•2  The  Coordinate  System 

We  have  to  determine  the  potential  functions  $  (  r ,  -  ,z  ) ,  r ,- ,z, ) 
and  ;(r,9,z^)  which  are  the  solution  of  Laplace's  equation  in  the  spaces 
(r,9,z^),  (r.d.z^)  and  (r,9,z^)  respectively.  Since  we  have  a  problem  in 
which  the  natural  boundaries  are  spheroids  with  center  at  the  origin  of 
coordinates,  it  is  suitable  to  employ  [5-9!  the  spheroidal  coordinate 
systems  defined  by 
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where  q j  end  p,  ere  parasmters  which  cen  be  determined  for  any  point  whose 
coordinates  (r,s)  ere  known.  The  ere  constants  to  be  determined  later. 

Let  Oj  denote  the  value  q^  on  the  spheroidal  surface,  then  che  chree 
coordinate  systems  coincide  on  the  surface  of  spheroid  If  the  following 
equalities  are  satisfied 
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It  is  interesting  to  note  that  when  q.  “  ■  Eq.  (4.2.2)  reduces  to 

the  equation  of  the  spheroidal  surface  described  by  Eq.  (4.1.1)  and  the 

2  z 

parameters  p.  become  independent  of  and  equal  The  coordinate  systems 
defined  by  Eq.  (4.2.2)  can  be  represented  graphically  in  the  physical  space 
for  real  values  of  vT.  For  example,  when  the  cavity  is  spherical  and  ,  = 

0.5  and  2.0,  Eq.  (4.2.2)  generates  two  types  of  spheroidal  and  hvperbo- 
loidal  surfaces  that  have  a  common  spherical  surface  corresponding  to  the 
same  value  ^ .  These  coordinate  systems  are  shown  in  Figs.  5  and  b. 

4 . J  The  Potential  Functions 


As  indicated  in  Section  3.2,  the  potential  functions  <t> ,  ,  and  y 

must  satisfy  Laplace's  equation  in  tne  spaces  (r,-),z,),  (r.-'.z,)  and 

(r,-,i  )  respectively.  It  is  well  known  [20,  29,  34 [  that  the  products 

P  m(p.)  Qm  (q.)  (cos  m^,  sin  m>'-j  are  harmonic  in  the  spaces  ( r ,  ■'■  ,  z  )  and 

n  )  n  )  : 

-rn  ^ 

regular  In  the  region  exterior  to  a  spheroid,  where  P  (d  )  and  <)' "(q  )  are 

n  i  n  ’ 

the  associated  Legendre  functions  of  the  first  and  second  kind  respectively 
and  will  be  defined  explicitly  later*  For  the  problem  at  hand,  the  poten¬ 
tial  functions  are  constructed  by  the  following  combination  of  the  Legendre 
associated  functions  of  first  and  second  kind: 
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because  they  are  harmonic  In  (r.e.Xj)  and  their  first  and  second  deriva¬ 
tives  vanish  at  Urge  distances  froe  the  spheroid.  Furthermore,  as  will  be 
shown  later,  the  first  derivative  of  the  chosen  potential  functions  with 
respect  to  r  or  x^  can  be  written  in  sleple  fores  which  lead  to  possible 
solutions  of  the  first  and  second  boundary  value  problems. 

In  these  potential  functions,  P^tpj )  and  C^Cq^  >  are  the  associated 
Legendre  functions  of  degree  n  and  order  a,  of  the  first  and  second  kind 
respectively.  They  are  defined  (20,27,34)  by 
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Legendre  associated  functions  defined  by  Eqs.  (4.3.2)  and  (4.3.6)  are 
listed  in  Appendix  A  for  different  values  of  n  and  ra. 

The  Legendre  associated  functions  defined  by  Eqs.  (4.3.2)  and  (4.3.7) 
satisfy  the  following  Identities 
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Th«  first  derivatives  of  Pn°(p)  and  Qn™(q)  with  respect  Co  che  argument  are 
given  by 
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(4.3.10) 


(q2  -  l)1  Q^nrt’1^  (q)  +  m  q  Q™(q) 


The  representation  of  the  potential  functions  in  terms  of  the  associ¬ 
ated  Legendre  function  as  shown  in  Eq.  (4.3.1)  enables  us  to  write  the 
partial  derivative  with  respect  to  r  and  in  very  simple  forms.  Recall¬ 
ing  Eq.  (4.2.2),  one  can  easily  obtain: 
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(4.3.11) 


®T  •*-*ns  of  Eqs.  (4. 3.8)-(4.3. 10) ,  it  can  be  readily  shown  that 


cos  m6  +  B  sin  m0 ] 


(4.3.12) 
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T~  *  l  l  2  p"m(P^  Q„Cq.)  (A,  c°s  m6 +B  sin  m9] 
3zj  n-0  m-0  n  J  n  j  jnm  jam 


(A. 3. 13) 


and  Eq.  (4.3.1)  can  also  be  written  in  the  form 


(4.3.14) 


[A  cos  ad  +  B,  sin  me] 
jnm  jnm 


4.4  Some  Remarks  on  the  Potential  Functions 


(i)  The  Legendre  associated  function  of  the  first  kind  used  in  the 
expressions  of  the  potential  functions  in  Eq.  (4.3.2)  is  equivalent 
to  the  following  definition 


where 


p>>  •  <-»■ » 


p“(p)  -  (-i)m  d-P2)m/2  ~ap-n-(P- 

n  dp 


(4.4.1) 


(4.4.2) 


for  all  values  of  m  <  n. 


(ii)  For  values  of  m  >  n,  Eq.  (4.3.2)  generates  a  function  of  p  singular 
at  p  *  1  or  p  =  -1.  This  situation  is  encountered  in  our  potential 
functions  when  m  ■  n  and  m  =  n+1.  However,  it  will  be  shown  later 
that  for  these  two  particular  cases,  singularities  inside  the  medium 
are  always  removable. 


(iii)  The  sum  with  respect  to  m  in  the  expression  of  the  potential  func¬ 
tion  is  truncated  after  the  value  m  *  n+1. 
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CHAPTER  5 


THE  FIRST  BOUNDARY  VALUE  PROBLEM 


5.1  Introduction 

The  first  boundary  value  problem  consists  of  finding  stresses  and  dis¬ 
placements  of  an  elastic  body  In  equilibrium  when  the  body  forces  are  known 
and  the  displacements  of  the  body  are  prescribed.  In  this  chapter,  the 
first  boundary  value  problem  Is  solved  for  an  unbounded  elastic  trans¬ 
versely  Isotropic  medium  containing  a  spheroidal  cavity  with  zero  body 
forces  and  vanishing  stress  and  displacement  fields  at  Infinity.  Explicit 
solutions  are  presented  for  two  problems  involving  a  rigid  spheroidal 
inclusion. 

5.2  Fundamental  Formulae 

The  components  of  the  displacement  vector  in  Eqs.  (3.2.13)  can  be 
written  as 


u 

r 

=  url 
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U  „  +  U  , 

r2  r3 
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3  U  , 
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u  «  +  u  ^ 

(5.2.1) 
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Substituting  the  expressions  of  and  ij>  defined  by  Eqs. 

Eqs.  (5.2.2)  we  obtain: 
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(m-1)  .  (m-1). 


[A.  cos  m6  +  B,  sin  me’ 
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(m-1) 


[B.  cos  m6  +  A,  sin  m6] 
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On  the  surface  of  a  spheroidal  cavity  (q^  -  and  -  p),  the 
of  the  displacement  vector  become: 


(5.2.2) 


.3.1)  into 


(5.2.3) 


components 
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2k.  2k,  m 

,  -  — -  Qm(c.)  A.  4  —  Q®(c,)  A 
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(5.2. 5) 
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(nrfl) , 
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xn  1  lnm  n  2  2nm  n  3 
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With  respect  to  rectangular  Cartisian  coordinates,  the  components  of 
the  displacement  vector  (ux.  uy,  uz )  can  be  written  in  terms  of  (ur,  u? , 
uz)  as  follows: 
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u  ■  u  cos  0  -  uAsin  6 
x  r  9 


u  ■  u  sin  0  +  u„cos  0 

y  r  0 


(5.2.6) 


u  *  u 
z  z 


Substituting  Eqs.  (5.2.4)  into  Eqs.  (5.2.6)  we  obtain: 
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(5.2.7) 
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If  a  displacement  vector  on  the  surface  of  the  spheriodal  cavity  is 
prescribed,  then  the  components  of  that  vector,  in  Cartesian  coordinates, 
u,  v,  and  w  can  be  represented  by  a  series  expansion  of  spherical  harmonics 
as  follows: 
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/  (u  cos  m2  +  u  sin  me)  P  (p) 
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(5.2.8) 
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where  P^Cp)  is  the  Legendre  associated  function,  defined  by  Eq.  (4.3.2), 
satisfying  the  following  orthogonality  condition 


in  which  '0  *  2  and  3m  -  1  for  m  4  0.  The  coefficients  vnm,  vnm,  wnm,  and 
wnn)  can  be  evaluated  in  a  similar  manner.  Equating  coefficients  in  the 
identical  trigonometric  and  spherical  functions  in  (5.2.7)  to  those  in 


(5.2.8)  yields 
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(5.2.11) 
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forms  as 
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In  which  [R^]  and  (Rnm]  are  3x1  aatrices  defined,  for  different  values 
of  m,  as  follows: 
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for  m  =  n  +  1 
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o 
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E**nn+1] 

1  (vnn+%n} 

L  o  j 

0  j 

(5.2. 13. d) 


If  m  ■  0,  the  stress  field  becomes  Independent  of  0  .  This  is  the 

axisymmetric  case  which  will  be  discussed  in  the  next  section.  Equations 
(5.2. 1 2a )  and  (5.2.12b)  have  a  unique  solution  if  the  determinant  A  is  not 
zero: 
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(5.2.14) 


v  n  1  n 


As  mentioned  in  the  previous  chapter,  the  Legendre  associated  function  of 
the  first  kind  defined  by  (4.3.2)  is  singular  at  p  =1  and  -1  for  the  two 
cases  corresponding  to  m  »  n  and  ra  *  n+I .  However,  the  equations  ann  = 

ann+l  “  ann  *  ann+l  *  Ynn+1  =*  Ynn+i  *  0  automatically  remove  these 
singularities  inside  the  medium. 


5.3  The  Axisymmetric  Problems 

Consider  the  class  of  problems  in  which  the  potentials  do  not  depend 
upon  6,  therefore  the  first  boundary  value  problem  can  be  deduced  directly 
from  the  preceeding  section  by  taking  ra  »  0  for  which  the  potential 


functions  become: 
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(5.3.1) 


(5.3.2) 


and  Che  displacements  in  cylindrical  coordinates  take  che  form: 


vr  j0  -2  v  p;;S>  «J(y 


’■)  A-in0  PnCpj 2  Q"(V 


“•j-° 


The  components  of  che  prescribed  displacement  vector  on  che  surface  of  the 
spheroidal  cavity,  in  che  axisymmetric  case,  can  be  represented  in  the 
forms: 


U  *  l  \l  Pn  (P} 
n*0 


”  *  Jo  W-  P”<P) 


(5.3.3) 


Here  u  and  w  are  the  components  of  the  displacement  vector  in  the  r-  and 
z-direction.  Expressing  Eqs.  (5.3.2)  on  the  surface  of  the  spheroid  and 
equating  their  coefficients  to  those  In  Eqs.  (5.3.3)  we  obtain: 
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(5.3.4) 


°.v  s-!  !v:""  we  determine  the  unknowns  A,  nf>  and  A->nn  which  enable  us 

to  determine  the  stress  field  without  difficulty. 

To  illustrate  the  method  developed  in  this  chapter,  we  will  consider 
the  problem  of  a  rigid  spheroidal  inclusion  embedded  in  a  transversely 
isotropic  medium  when  the  medium  is  subjected  to  the  following  loads  at 
large  distances  from  the  inclusion: 

(I)  Uniaxial  tension  in  the  direction  of  the  axis  of  symmetry  of  the 
medium. 

(II)  Pure  shearing  stress  in  the  plane  perpendicular  to  the  axis  of 
symmetry  of  the  medium. 


5 . 4  Rigid  Spheroidal  Inclusion  Under  Axlsymmetrlc  Uniaxial  Tension 
The  boundary  conditions  in  this  case  are: 
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(5.4.1) 


u  ■  u 
r 


u 

z 


0  on  the  surface  of  spheroid 


(5.4 


2) 


In  the  absence  of  the  inclusion,  the  uniform  stress  field  (Eqs.  5.4.1)  can 
be  extended  throughout  the  space,  thereby  violating  the  boundary  conditions 
in  Eqs.  (5.4.2).  The  stress  and  displacement  field  are  found  to  be: 
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We  now  see*  a  solution  union,  upon  superposition  on  the  umiorm  iLeld, 
removes  the  residual  displacements,  given  by  Eqs.  (3.-*.-*),  on  the  sur race 
of  the  spheroid.  Therefore,  consider  the  problem  where  the  negative  of  the 
displacements  given  by  tqs.  (5. a. 4;  are  specified  on  the  boundurv  of  the 
spheroidal  inclusion.  These  displacements  are: 
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in  which  a  and  b  are  the  semiaxes  of  the  spheroid.  From  Eqs.  (3.3.3)  it  is 
easily  found  that  the  only  non-zero  coefficients  are  those  corresponding  to 


n  =  1.  Thus, 
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Equations  (5*3.4)  become  in  this  case: 
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Solving  the  simultaneous  equations  we  obtain: 
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The  potential  functions,  therefore,  are 
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Substituting  Eqs.  (5.4.9)  Into  Eqs.  (3.2.13)  and  (3.2.14;  and  adding  the 
solution  obtained  in  the  absence  of  the  inclusion,  we  obtain  the  desired 


solution  for  which  the  displacements  and  stresses  become: 
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5.5  Rigid  Spheroidal  Inclusion  Under  Pure  Shearing  Stress  in  the  Plane  of 


Isotropy 


The  boundary  conditions  in  this  case  are: 
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(5.5.1) 


ux  ■  Uy  ■  uz  «  0  on  the  surface  of  the  spheroid 


(5.5.2) 


the  loading  conditions  (5.5.1)  is  equivalent,  in  cylindrical  coordinates. 
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In  the  absence  of  the  inclusion,  Eqs.  (5.5.1)  or  (5.5.3)  represent  the 
stress  field  in  the  medium.  The  displacement  field,  accordingly,  is  found 
to  be : 
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Consider  the  problem  where  the  negative  of  these  displacements  are 
specified  on  the  boundary  of  the  spheroid.  Thus, 
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Comparing  Eqs.  (5.5.5)  to  Eqs.  (5.2.6),  we  find  that  the  only  non-zero 
coefficients  are: 
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By  inspecting  Eqs.  (5.2.11)  and  (5.2.13)  we  conclude  that  the  desired 


problem  is  corresponding  to  n  =1,  m  *  2,  A .  *  0,  and  B.  1 0  to  be 


determined.  Therefore,  Eqs.  (5.2. 1 2 B )  become: 
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which  have  the  solution: 
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The  potential  functions  corresponding  to  the  problem  at  hand  are: 
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Substituting  into  Eqs.  (3.2.2)  and  adding  the  solution  of  that  in  the 
absence  of  the  cavity,  we  obtain  the  following  displacements  and  stresses: 
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At  the  spheroid  interface,  these  equations  may  be  reduced  to: 
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If  we  obtain  a  transformation  of  the  stress  components  from  cylindrical 
coordinates  to  spheroidal  coordinates  by  means  of  the  equations  in  Appendix 
B,  we  obtain 
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where  D  is  defined  by  Eq.  (5.5.8)  and  s  •  — . 


In  a  special  case  where  the  medium  is  isotropic,  Eqs.  (5.5.12)  become 
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Kir  i  r  1  i ti  spherical  inclusion  we  must  evaluate  the  limit  of  Dq  as 
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Substituting  the  limiting  value  of  Dq  into  Eqs.  (5*5.13),  we  obtain  the 


CHAPTER  6 


THE  SECOND  BOUNDARY  VALUE  PROBLEM 


6.1  Introduce  ion 

The  second  boundary  value  problem  consists  of  finding  stresses  and 
displacements  of  an  elastic  body  in  equilibrium  when  the  body  forces  are 
known  and  the  surface  forces  are  prescribed.  In  succeeding  sections  of 
this  chapter  the  second  boundary  value  problem  is  solved  for  equilibrium  of 
an  elastic  transversely  Isotropic  medium  containing  a  spheroidal  cavity. 
Examples  are  given  for  a  variety  of  constant  loadings  applied  at  suffi¬ 
ciently  large  distances  from  the  unloaded  cavity  In  the  absence  of  body 
forces . 


6.2  Fundamental  Formulae 


Let  nr  and  nz  denote  the  components  of 
cosines)  to  the  spheroidal  surface  for  which 
evaluated  from  the  expressions 
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The  projections  tr,  t2,  and  C-  of  che  traction  vector  acting  upon  an 
area  to  which  the  normal  is  given  by  the  direction  cosines  nr  and  nz  can  be 
expressed  as  follows 
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We  shall  express  these  components  in  terms  of  the  potential  functions 
defined  by  Eqs.  (4.3.1).  First,  however,  we  will  obtain  an  expression  that 
will  be  used  frequently  in  the  analysis: 
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On  the  surface  of  the  spheroidal  cavity  we  have  q^  »  Cj  and  p.  »  p.  Thus, 
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Using  Eqs.  (4.2.3)  and  (6.2.1)  we  obtain 
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It  is  noc  difficult  to  show  that  by  means  of  Eqs.  (3.2.8),  (o.2.o),  and 

(6.3.1),  Eqs.  (o.3.2)  cake  cne  form: 
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6.4  The  Components  of  the  Traction  Vector  Due  to 


The  components  of  Che  displacement  vector  due  to  the  potential  func- 
don  ij)  (the  "curl'*  portion)  can  be  written  as: 
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It  can  be  shown  that  Eqs.  (6.4.1)  satisfy  the  following  equation: 
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Subscituting  Eqs.  (3.2.1)  Into  Eqs.  (6.2.2)  and  setting  uZ3  ■  G,  Che  coopo- 
nencs  of  Che  traction  vector  due  to  the  potential  function  i|i  become: 
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By  using  Eqs.  (3.2.8),  (6.2.6),  (6.4.1),  and  (6.4.2),  it  can  be  shown  that 
Eqs.  (6.4.3)  have  the  form: 
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6 .  5  The  Components  of  the  Traction  Vector  t  •>  Terms  of  the  Legendre 
A b s •  ^ciated  Functions 

Expressions  of  che  displacement  components  are  given  by  Eqs.  (5.2.3) 
Substituting  these  expressions  into  Eqs.  (o.3.3)  and  (6. 4.1)  we  obtain 
after  some  algebraic  manipulation,  che  following  expressions  for  the  compo 
neats  of  the  traction  vector  on  the  surface  of  the  spheroidal  cavity: 
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By  setting 
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and  using  tha  following  transformation: 
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If  a  surface  force  on  Che  spheroidal  cavlcy  is  prescribed,  Chen 
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nents  of  chat  veccor,  in  Carcesian  coordinaCes,  t* ,  c^ ,  and 
represenced  by  a  series  expansion  of  spherical  harmonics  as 


(6.5.5) 


(6.5.6) 


Che  compo- 

0 

Cz  can  be 
follows : 


60 


oc  n 

01-“/  _  -  .  -m,  . 

t  *  t  1  (g  cos  mo  +  g  sm  mo)  P  (p) 

x  X  -  -  °nm  ®nm  n 

n«0  m»0 


0  1  -  - 

t  *  —  '  (h  cos  ra6  +  h  sin  m0)  P  m(p) 

y  '  -n  nm  nm  n 

n-0  ra-0 


0  1  —  ^ 
c  *  —  '  '  (l  cos  mo  +•  i  sin  0)  P  (p) 

z  ■  .  nm  nm  m  n 

n-0  m-0 


(6.5.7) 


in  which  the  coefficients  Snm*  §nm*  ^nm*  h^m,  inm*  and  inm  can  be  evaluated 
by  using  the  orthogonality  property  which  leads  to  equations  similar  to 
chose  of  Eqs.  (5.2.10).  By  equating  the  coefficients  in  the  Identical 
trigonometric  and  spherical  functions  in  (6.5.5)  to  those  in  (6.5.7)  we 
obtain  the  following  equations: 
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where  [vntBj  and  [  vnmj  are  defined  as  follows: 
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For  values  of  at  »  n  and  m  *  a  +  1 ,  singularities  at  p  •  1,  -1  are  removed 
by  the  equations  8nm  -  3nm  -  ynn+1  -  ynnn  -  ;nn+1  -  Cnn+i  "  0. 

6.6  The  Axisvmmetric  Problem 


For  the  axisvmmetric  problem  In  which  the  stress  field  is  independent 
of  6,  the  potential  functions  take  the  form 
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If  a  surface  force  on  the  spheroidal  cavity  is  prescribed,  the  compo¬ 
nents  of  this  vector  can  be  represented  by: 
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displacement  field  can  be  readiiv 
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As  an  Illustration,  we  will  consider  in  the  next  sections  the  problem 
of  transversely  Isotropic  medium  containing  a  spheroidal  cavity  when  the 
medium  is  subjected,  at  large  distances  from  the  cavity,  to 

(i)  Uniaxial  tension  in  the  direction  of  axis  of  symmetry  of  the 
medium  (z-di rect ion). 

(ii)  Hydrostatic  tension  in  the  plane  of  isotropy  (xy-plane). 

(iii)  Pure  shear  stress  in  the  plane  of  isotropy  (xy-plane). 
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nents  are  specified  on  che  spheroidal  cavitv,  therefore: 
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Comparing  Eqs.  (6.7.5)  to  Eqs.  (6.6.2),  we  find  that  the  only  oon-cero 
coefficient  is  l^g  *  -T0b.  Therefore, 
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and  Eqs.  (6.6.3)  which  determine  A^q  anc*  ing  become 
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from  which  we  obtain 
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The  pocenclal  functions  in  this  case  are: 
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Substituting  (6.7. 10)  into  (3.2.1)  and  adding  the  solution  of  the  problem 
in  the  absence  of  the  cavity,  the  displacement  and  stress  field  for  the 
problem  at  hand  become: 
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Ac  Che  spheroid  Incerface,  che  equations  of  che  stresses  In  the  medium  take 
Che  form: 
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By  transforming  Eqs.  (6.7.12)  from  cylindrical  coordinates  to  spheroidal 
coordinates  by  means  of  Eqs.  (B.3)  in  Appendix  B  we  obtain: 
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For  an  isotropic  medium  containing  a  spheroidal  cavity,  Eqs.  (6.7.13)  are 
reduced  by  the  limiting  process  to: 
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where  y  is  the  meridional  ant  it  shown  in  Fig.  8. 
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6.8  Hydrostatic  Tension  in  t!ie  Hlane  of  Isotrot 


The  boundary  conditions  of  this  problem  are 
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In  the  absence  of  the  cavity,  the  state  of  stress  is  represented  by  Eqs. 
(6.8.1).  Therefore,  the  tractions  that  must  be  removed  to  satisfy  the 
boundary  conditions  of  (p.8.2)  are: 
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We  now  consider  the  problem  in  which  the  negative  of  these  tractions  ire 
applied  on  the  surface  of  the  spheroidal  cavity.  Thus, 
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It  follows  that  the  only  non-zero  term  in  Eqs.  (6.6.2)  is  *  -2T0a. 
This  tera  is  corresponding  to  n  *  1.  The  coefficients  and  A^iq  can  b* 
determined  by  solving  the  equations 
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where  Hi  ,  Hg^g  ,  S110  •  and  5->  i  g 
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The  potential  functions  for  this  case  take  the  form: 
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The  stress  field  becomes: 
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On  the  spheroidal  surface,  these  equations  take  the  fora: 
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By  transforming  Eqs.  (0.8. 9)  from  cylindrical  coordinates  to  spheroidal 
coordinates  by  means  of  Eqs.  (3.3)  in  Appendix  3  we  ootam: 
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For  Isotropy,  by  caking  the  proper  Limit  we  find 
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For  an  isotropic  medium  containing  a  spherical  cavity,  Eqs.  (6.8.13) 
become 
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where  v  is  the  meridional  angle  shown  in  Fig.  8. 


6.9  Pure  Shear  Stress  in  the  Plane  of  Isotroc 


(6.8.14) 


The  boundary  conditions  are 
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In  the  absence  of  the  cavity,  the  state  of  stress  in  the  iaedium  can  be 


represented  by: 


(6.9. 3 ) 


or,  equivalently  in  cylindrical  coordinates: 
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therefore,  the  problem  at  hand  is  corresponding  to  the  case  in  which 
n  =•  1,  and  A  =*  0.  The  coefficients  B  are  determined  from: 
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The  solution  of  Eqs.  (6.9.9),  after  some  algebric  manipulation,  is  given  by 
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The  potential  functions  for  this  case  arc: 
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Substituting  Eqs.  (6.9.13)  into  Eqs.  (3.2.1)  and  adding  the  solution  of  the 
problem  in  the  absence  of  the  cavity  we  find, 
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On  the  spheroidal  surface,  the  stress  components  in  the  medium  take  the 
form: 
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The  stress  components  In  spheroidal  coordinates  are  obtained  by  using  the 
transformation  In  Appendix  B  and  take  the  form: 
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For  an  Isotropic  medium,  Eqs.  (6.9.16)  aay  be  reduced  to: 
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Furthermore,  for  an  isotropic  medium  containing  a  spherical  cavity  Eqs. 
(6.9.17)  are  reduced  to: 
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6.10  Uniaxial  Tension  in  the  x-Dlrection 


The  distribution  of  the  stress  field  throughout  a  medium  subjected  to 
uniaxial  tension  in  the  x-direction  T0  may  be  achieved  by  superposition  of 
the  solutions  for  the  following  two  cases: 

(i)  ax  -  Oy  ■  "2  T0 
(ii)  m  ~°y  m  ~2  To 

The  solution  of  case  (i)  may  be  obtained  directly  from  Section  (6.8) 
after  replacing  T0  by  T0.  On  the  other  hand,  the  solution  of  case  (ii) 
can  be  obtained  from  the  solution  of  Section  (6.9)  after  replacing  6  by 
(0  +  ~  '  and  x0  by  tq.  Therefore,  the  stress  components  at  the  spher¬ 
oidal  cavity  surface  for  the  desired  problem  are: 


v3  (l+k1) (l+kj)  4  2 


2sFl  (k2-kx) 


(r  +  s  **■)  ( 


/1  2  2i 

(1  -  s  v^) 


v^(r2  + s4  z2)Q^(o1)  F3 


J  (-1  + 


(l+k2)(l-v‘) 


,,  ,  ,  2  2  4  2 

(k1-k2)  r  +vLs  z 


]  cos  26 


2  2 

(1  -  s  v2) 

,  2  2  4  2.1,  7 

v2  r  +  v2  s  2  ^1  p2 

(l+k^d-v2)  r2 

~  “  2  2  4  2 

(k1~k2)  r  +v2s  z 
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!  CHAPTER  7 

t 

t 

I 

\  THE  EFFECT  OF  ANISOTROPY  ON  THE  STRESS  CONCENTRATION  FACTORS 

|  7,1  Introduction 

In  this  chapter,  numerical  results  will  be  presented  for  a  transverse¬ 
ly  isotropic  medium  containing  a  spherical  cavity  subjected  to  the 
following  load  cases,  applied  at  large  distances  from  the  cavity. 

(i)  Uniaxial  tension  in  the  direction  of  the  axis  of  symmetry  of  the 
material  (z-direct ion) . 

(ii)  Uniaxial  tension  in  the  direction  perpendicular  to  the  axis  of 
symmetry  of  the  material  ( x-direct ion) . 

Elementary  dimensional  considerations  show  that  the  stress  concent  ration 
factor  must  depend  upon  four  dimensionless  ratios  of  either  the  five  elas¬ 
tic  constants  (stiffnesses),  or  the  compliances,  or  the  engineering 
constants.  Adopting  the  engineering  constants  and  noting  that  Poisson's 

ratios  v  and  v  are  dimensionless,  the  remaining  two  dimensionless  ratios 

-  EC 

among  E,  E,  and  G  can  be  taken  as  —  and  ^  although  other  ratios  are 

possible.  To  show  the  effect  of  anisotropy  on  the  stress  concentration 
factor  in  the  vicinity  of  a  spherical  cavity,  we  consider  hypothetical 
materials  for  which  one  of  the  four  chosen  ratios  can  vary  while  the 
remaining  three  ratios  are  kept  fixed  at  values  corresponding  to  an  isotro¬ 
pic  material  whose  Poisson's  ratio  equals  0.25.  Therefore,  the  ratios 
have,  at  isotropy  with  v  =  0.25,  the  following  values 

v  =  v  =  0.25  ,  |  =  1.0  ,  |  = 


0.4 


(7.1.1) 


Throughout  this  chapter,  we  define  the  stress  concentration  factor  as  the 
ratio  of  the  maximum  principal  stress  on  the  surface  of  the  cavity  to  the 
magnitude  of  the  uniform  stress  field  applied  at  infinity.  Positive  and 
negative  signs  are  used  to  denote  tensile  and  compressive  stresses 
respect i vely . 

Numerical  results  for  the  stress  concentration  factor  are  plotted  on  a 
coordinate  system  in  which  the  abscissa  represents  the  varying  ratio  and 
the  ordinate  represents  the  stress  concentration  factor  denoted  by  Kc  and 
Kc  for  tensio;  and  compression  respectively.  Negative  values  of  Poisson's 
ratios  do  theoretically  exist  but  they  may  not  be  physically  attainable. 
Therefore,  they  are  omitted  from  the  plots  and  discussion.  The  asterisk 
(*)  on  the  plots  indicates  the  upper  or  lower  limit  of  any  of  the  four 
ratios  for  positive  definiteness  of  the  strain  energy  function  (Eq.  2.4.3) 
when  the  other  three  ratios  are  fixed.  For  each  value  of  the  varying 
ratio,  the  locations  of  Kt  and  Kc  are  calculated  and  shown  on  a  plot  in 
terms  of  two  angles  y  and  0  shown  in  Fig.  8. 


7.2  Uniaxial  Tension  in  the  z-Direction 


This  problem  was  first  investigated  by  Chen  [5]  who  evaluated  the 
stress  concentration  factor  for  a  few  transversely  isotropic  materials.  In 
this  section,  we  will  present  more  numerical  data  than  has  been  reported 
previously  in  the  literature. 

It  is  well  known  that  the  tensile  and  compressive  stress  concentration 
factors  for  a  spherical  cavity  in  an  isotropic  medium  equal  and 
-  y  respectively.  The  former  occurs  on  the  equatorial  line  and 


changes  slightly  from  1.929  to  2.167  as  Poisson's  ratio  traverses  the 
range  0  to  0.5  and  the  latter  occurs  at  the  pole  and  changes  from  -0.214  to 
-1.1  lo  as  .  traverses  the  range  ;)  to  0.5. 

We  have  calculated  the  tensile  and  compressive  stress  concent r at  ion 
factors  for  a  wide  range  of  crystalline  and  non-crystalline  transversely 
isotropic  materials  Including  some  composite  materials  possessing  high 
Young's  modulus  in  the  direction  of  the  fibers.  Results  of  the  computa¬ 
tions  are  presented  in  Table  5  which  shows  a  substantial  increase  in  the 
maximum  tensile  stress  for  highly  anisotropic  materials  (i.e.,  Graphite 
Thornel).  Furthermore,  the  highest  tensile  and  compressive  stresses  for 
all  of  the  anisotropic  materials,  listed  in  Table  5,  occur  at  the  equatori¬ 
al  line  and  pole  respectively.  This  is  a  well  recognized  situation  in 
isotropic  materials  but  this  observation  should  not  lead  to  a  general  con¬ 
clusion  applicable  to  other  anisotropic  materials.  To  provide  further 
insight  into  the  influence  of  anisotropy  on  the  stress  concentration 
fictor,  sensitivity  analyses  for  the  four  ratios  mentioned  in  the  proceed¬ 
ing  section  have  been  made  and  will  be  presented  next. 

7.2.1  Effect  of  — 

Figure  9  shows  the  variation  of  the  nighest  tensile  and  compressive 

E 

stress  factors  Kt  and  Kc  against  the  variation  of  the  ratio  — .  It  is  seen 

1  E 

that  for  values  of  —  <  —  1,  Kt  changes  slightly  and  remains  close  to  2.0 

6  E 

E  EE 

but  for  —  >  1,  Kt  increases  appreciably  as  —  increases  (i.e.,  for  -  =  10 

E  EE 

E 

and  100,  Kc  »  4.214  and  ll.ll).  Conversely,  the  factor  Kc  decreases  as  ^ 

E 

increases  and  becomes  negligible  for  very  large  values  of  (i.e.,  —  =  10 

and  LOO,  Kc  =  -0.140  and  -0.041). 


'  n 


i  ■  t-M  «.>  |.| 


In  Fig.  10,  ic  is  shown  that  for  —  =  0.167  the  highest  tensile  stress 

E 

is  calculated  at  y  =*  49.3°  while  the  highest  compressive  stress  is  calcula- 

a  E 

ted  at  v  »  19. 2 3 .  As  —  deviates  and  becomes  larger  than  0.167,  the  loca- 

E 

tion  of  the  highest  tensile  stress  starts  shifting  toward  the  equatorial 

E 

line  (Y  =  90°)  and  remains  ther  for  —  >  0.5,  whereas  the  location  of 
highest  compressive  stress  starts  shifting  toward  the  pole  and  remains  at 

E 

the  pole  for  values  of  p  >  0.3. 

E 


7.2.2  Effect  of 


The  variation  of  the  tensile  and  compressive  stress  concentration 

Q 

factors  with  —  are  shown  in  Fig.  11.  It  can  be  observed  that 
E 

Q 

(i)  As  the  ratio  —  rises  from  0.01  to  1,  Kt  decreases  from  7.400  to 
1.749. 

G  G 

(ii)  Neither  Kt ,  for  —  >  1,  nor  Kc,  for  —  >  0.01,  seems  to  be  sensi- 

G 

tive  to  the  variation  of  — . 

E 

Q 

(iii)  For  a  sufficiently  large  value  of  — ,  Kt  and  Kc  approach  the 
values  1.913  and  -0.633  respectively. 

Figure  12  shows  that  the  location  of  the  highest  tensile  stress  is 

Q 

found  to  be  on  the  equatorial  line  (y  =  90°)  for  values  of  —  _<  0.8.  This 

location  starts  shifting  toward  the  pole  as  increases  and  deviates  from 
r 

0.8.  When  —  becomes  sufficiently  large,  the  highest  tensile  stress  occurs, 
E 

approximately,  at  y  -  57.3°.  It  should  also  be  noted  th3t  for  all  vilif- 

r 

of  — ,  the  highest  compressive  stress  occurs  at  the  pole. 
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7.2.3  Effect  of  Poisson's  Ratios  v  and  v 


Figures  13  and  14  show  clearly  that  Poisson's  ratios  v  and  v  have 
little  effect  on  the  stress  concentration  factors  Kt  and  Kc.  When  v  trans- 
verses  the  range  0  <  v  <  J—  »  0.612,  Kt  changes  slightly  between  2.028  and 
2.160  and  Kc  Increases  from  -0.333  to  -0.903.  For  values  of  0  <  -  <  ^ 
plotted  in  Fig.  14,  Kt  appears  to  be  insensitive  to  the  variation  of  v  and 
remains  close  to  2.0  whereas  Kc  increases  from  -0.45  to  -2.80.  The  loca¬ 
tions  of  the  highest  tensile  and  compressive  stress  are  found  to  be  at  the 
equator  and  pole  respectively  for  all  values  of  v  and  v. 


7.3  Uniaxial  Tension  in  the  Direction  Perpendicular  to  the  Axl9  of  Elastic 
Symmetry  (x-Direction) 


Obviously,  for  the  case  of  an  Isotropic  body  with  a  spherical  cavity, 
the  solutions  corresponding  to  tensions  in  the  x-and  z-direction  are  equi¬ 
valent  since  these  two  directions  are  elastically  and  geometrically 
identical. 

For  transversely  isotropic  materials,  the  magnitude  of  the  highest 
principal  tensile  and  compressive  stresses,  Kt  and  Kc,  are  calculated  for  a 
variety  of  materials.  Results  of  computation  are  given  in  Table  6.  It  is 
found  that  for  the  majority  of  these  materials,  the  highest  principal 
tensile  stress  occurs  at  the  pole  (  y  -  0°).  However,  for  some  materials 
indicated  by  an  asterisk  in  Table  6,  the  highest  principal  tensile  stress 
occurs  away  from  the  pole.  On  the  other  hand,  the  highest  principal 
compressive  stress  for  all  the  materials  in  Table  6  occurs  at  0  ■  0°  and 
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y  ■  90°.  The  effect  of  anisotropy  on  on  the  stress  concentration  factor  on 
the  stress  concentration  factor  of  the  problem  at  hand  is  discussed  next. 


7.3.1  Effect  of  — 
E 


Figure  15  presents  the  variation  of  the  maximum  principal  tensile  and 
compressive  stress  factors,  Kt  and  Kc,  on  the  cavity  surface  against  the 

E 

variation  of  the  ratio  — .  It  is  seen  that  the  factor  Kr  is  insensitive  to 

E  u 

E 

the  variation  of  — ,  changes  from  a  value  of  slightly  above  2.5  to  slightly 
E  __ 

£ 

below  2.0,  and  approaches  a  value  of  1.946  for  fairly  large  values  of  |r. 

E 

In  contrast,  the  factor  Kc  is  considerably  affected  by  the  variation  of 

E 

attains  its  absolute  minimum  (Kc  *  2.022)  at  —  ■  1,  and  increases  monotoni- 

E 

E 

cally  as  —  increases.  The  variation  of  the  location  of  the  maximum 

E 

principal  tensile  stress  when  —  varies  is  shown  in  Fig.  16  while  the 

E 

maximum  principal  compressive  stress  occurs  at  6  »  0°  and  y  *  90°. 


7.3.2  Effect  of  | 


l 


Q 

In  Fig.  17,  Kt  is  plotted  as  a  function  of  the  ratio  — .  It  is  seen 

C  G 

that  if  —  rises  from  0.01  to  0.4,  Kt  decreases  from  2.923  to  2.022.  When  ^ 

changes  from  0.4  to  9,  Kt  changes  slightly  from  2.022  to  1.893.  For  values 

G  G 

of  —  greater  than  9,  as  —  increases,  the  magnitude  of  Kt  increases.  The 

E  E 

location  at  which  Kt  is  evaluated  is  shown  in  the  same  plot.  A  similar 

plot  is  made  for  Kc  and  shown  in  Fig.  18.  It  appears  that  Kc  is  insensi- 

G  G  G 

tive  to  the  variation  of  —  when  <  10  but  it  increases  as  the  ratio 

EE  t 

increases  (i.e.,  ^  ■  100,  Kt  -  5.003  at  0  *  50°  and  y  ■  86°,  and  Kc  * 

-4.726  at  J  “  40°  and  y  *  86°). 
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7.3.3  Effect  of  v  and  v 

Fig.  19  shows  that  neither  Kt  nor  Kc  is  sensitive  to  the  variation  of 
v.  The  location  at  which  Kt  is  evaluated  is  shown  in  Fig.  20  whereas  Kc  is 
found  to  be  at  6  *  0°  and  y  “  90°  for  all  plotted  values  of  v.  In  Fig.  21 
Kt  and  Kc  are  plotted  as  a  function  of  v.  It  is  seen  that  when  v  changes 
from  0  to  0.875,  Kt  rises  from  2.002  to  4.017.  The  location  at  which  Kt  is 
evaluated  is  9  *  90®  and  y  »  90®  for  v  _<  0.25  and  is  6  *  90®  and  Y  *  0®  for 
v  _>  0.25.  In  the  same  plot,  Kc  changes  from  -0.525  to  —1.115,  when  v 
changes  from  0  to  0.875,  and  occur  at  6  ■  0®  and  Y  *  90®. 
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Table  5.  Stress  Concentration  Factor  on  the  Surface  of  Spherical  Cavity 
Under  Uniaxial  Tension  in  the  Direction  of  the  Axis  of  Elastic 
Symmetry 


Material 

Stress 

Concentration  Factor 

Tension  (Kt) 

Compression  (Kc) 

Beryllium 

1.960 

-0.248 

Bone  (fresh  phalanx) 

2.328 

-0.654 

Cadm iura 

1.76b 

-1.033 

Ceramics  (BaTiO^) 

2.048 

-0.729 

Cobalt 

2.285 

-0.690 

Ecologite 

2.130 

-0.532  ! 

Graphite 

2.753 

-0.669 

Graphite  Thornel  50 

4.683 

-0.146  j 

Graphite  Thornel  75 

5.637 

-0.H7 

J 

Hafnium 

2.1U7 

-0.616  | 

Ice 

2.231 

-0.664 

*. 

Magnes ium 

2.129 

-0.b33 

> 

Micha  Schist 

1.980 

-0.538  | 

Quartz 

2.025 

-0.491  v 

Rhenium 

2.220 

-0.596  !;• 

Silver  Aluminum 

2.210 

-0.743  jj 

Titanium 

2.155 

-0.704 

Zinc 

1.621 

-1.000  ij 

Zinc  Oxide 

2.19b 

-0.769  jj 

Zirconium 

2.243 

-0.677  ;S 

Isotropic  medium 
(  v  -  0.25) 

2.022 

-0.587  ^ 

8 

51 


Table  6.  Stress  Concentration  Factor  on  the  Surface  of  Spherical  Cavity 
Under  Uniaxial  Tension  Perpendicular  to  the  Axis  of  Elastic 
Symmetry 


Material 

Stress 

Concentration  Factor 

Tension  (Kc) 

Compression  (Kc) 

Beryllium 

1.905 

-0.278 

Bone  (fresh  phalanx) 

2.127 

-0.976 

Cadmium 

2.316* 

-0.763 

Ceramics  (BaTiO^) 

2.069* 

-0.7-24 

Cobalt 

2.262 

-0.770 

Ecologite 

2.476* 

-0.630 

Graphite 

2.918* 

-0.946 

Graphite  Thornel  50 

1.935 

-2.0b6 

Graphite  Thornel  75 

1.921 

-2.793 

Hafnium 

2.075 

-0.649 

Ice 

2.142 

-0.726 

Magnesium 

2.105 

-0.665 

Micha  Schist 

2.074* 

-0.518 

Quartz 

2.021* 

-0.526 

Rhenium 

2.146 

-0.566 

Silver  Aluminum 

2.176 

-0.812 

Titanium 

2.177 

-0.736 

Zinc 

2.315 

-0.637 

Zinc  Oxide 

2.186 

-0.792 

Zirconium 

2. 1 50 

-0.745 

Isotropic  medium 
( v  -  0.25) 


2.022 


-0.587 


CHAPTER  8 


j 

SUMMARY,  CONCLUSION,  AND  RECOMMENDATION  FOR  FURTHER  STUDY 

I 

I 

i  8.1  Summary  and  Conclusion 

The  principal  result  of  this  study  has  been  the  development  of  expli¬ 
cit  analytical  solutions  for  the  (non-axisymmetrlc)  first  and  second 
boundary  value  problems  of  elasticity  theory  for  a  spheroidal  cavity 
embedded  in  a  transversely  isotropic  medium.  The  analysis  is  based  upon 
solutions  of  the  homogeneous  displacement  equations  of  equilibrium  In  terms 
of  three  quasi-harmonic  potential  functions  taken  in  a  special  combination 
of  the  associated  Legendre  functions  of  the  first  and  second  kind.  Exact 
solutions  have  been  obtained  for  problems  involving  a  region  containing  a 
rigid  spheroidal  inclusion  and  a  region  containing  a  traction-free  cavity 
subjected  to  constant  loadings  applied  at  sufficiently  long  distances  from 
the  cavity.  A  tractable  problem  which  can  be  treated  in  a  similar  fashion 
is  the  hyperbololdal  notch  in  a  transversely  isotropic  material  under 
arbitrary  loadings  (see  Reference  [8]). 

8.2  Recommendation  for  Further  Study 

The  applications  of  the  present  approach  are  by  no  means  exhausted  in 
this  work.  Further  study  should  be  made  of  the  following: 


Vv-  -y1 
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1.  Problems  which  can  be  solved  directly  by  the  present  approach.  These 
Include  solutions  of  spheroidal  cavities  and  inclusions  under  loading 
conditions  different  from  those  we  have  obtained,  as  well  as  numerical 
investigations  including  the  effect  on  the  stress  concentration  factors 
of  the  shape  ratio  of  the  spheroidal  cavity  or  inclusion  and  the  effect 
of  anisotropy  on  the  decay  of  stresses  with  distance  from  the  cavity  or 
inclusion. 

2.  The  mixed  and  mixed-mixed  boundary  value  problems  for  a  transversely 
isotropic  medium  containing  a  spheroidal  cavity  or  inclusion. 


3.  Problems  in  which  the  elastic  medium  is  bounded  by  other  geometries 
(i.e.,  the  hyperboloidal  notch). 
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Figure  1.  Hexagonal  Arrays 


Figure  2.  Honeycomb  Structures 


Figure  3.  Stratified  Materials 
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figure  6. 
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Figure  9.  Variation  of  Stress  Concentration  with  Modulus  Ratios 
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Figure  10.  Meridional  Position  of  Stress  Concentration  versus 
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Figure  21.  Variation  of  Stress  Concentration  with  v 
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APPENDIX  B 


TRANSFORMATION  OF  THE  STRESS  COMPONENTS  FROM 
CYLINDRICAL  TO  SPHEROIDAL  COORDINATES 


Let  Che  spheroidal  coordinate  system  be  defined  by  means  of  the  trans¬ 


formation 


r  -  a(q2  -  l7(l  -  p2)'^ 


2  -  a  q  p 


(B.l) 


in  which  it  is  customary  to  use  q  ■  cosh  n  and  p  -  cos  0  where  the  ranges 
of  ",  t,  and  6  are 


0  <  n  <  ®  ,  0  <  #  <  it  ,  0<9<2tt 


(B  .2) 


and  let  ~q  be  the  value  on  the  surface  of  a  spheroid  whose  semi-axes  are 
a  and  b.  Then,  a  •  a  cosh  ng ,  b  ■  a  sinh  Hg ,  and  the  components  of  the 
stress  tensor  along  three  orthogonal  directions  (n,  $  ,  0 )  can  be  written  in 
terms  of  the  components  of  the  stress  tensor  in  the  cylindrical  coordinate 
system  as  follows. 


Inn  “  7777  r  +  3  2  °«  +  2s  rz 


1  .42  .2  .  2 

7777  3  2  °-  +  r  °~'29  r20~] 


°00  “  a99 


(B.  3) 
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